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Abstract 

A distributed detection problem over fading Gaussian multiple-access channels is considered. Sensors observe 
a phenomenon and transmit their observations to a fusion center using the amplify and forward scheme. The fusion 
center has multiple antennas with different channel models considered between the sensors and the fusion center, 
and different cases of channel state information are assumed at the sensors. The performance is evaluated in terms of 
the error exponent for each of these cases, where the effect of multiple antennas at the fusion center is studied. It is 
shown that for zero-mean channels between the sensors and the fusion center when there is no channel information 
at the sensors, arbitrarily large gains in the error exponent can be obtained with sufficient increase in the number 
of antennas at the fusion center. In stark contrast, when there is channel information at the sensors, the gain in 
error exponent due to having multiple antennas at the fusion center is shown to be no more than a factor of S/tt 
for Rayleigh fading channels between the sensors and the fusion center, independent of the number of antennas at 
the fusion center, or correlation among noise samples across sensors. Scaling laws for such gains are also provided 
when both sensors and antennas are increased simultaneously. Simple practical schemes and a numerical method 
using semidefinite relaxation techniques are presented that utilize the limited possible gains available. Simulations 
are used to establish the accuracy of the results. 

I. Introduction 

Sensors are becoming commonplace in factories, environmental, and home appliance monitoring, as 
well as in scientific study. In many such applications, a number of independent sensors each make a local 
observation, which are transmitted to a fusion center (FC) after limited initial processing at the sensors, 
and combined at the FC to calculate a global decision [IJ. Sensors may adopt either a digital or an 
analog method for relaying the sensed information to the FC. The digital method consists of quantizing 
the sensed data and transmitting the digital data over a rate-constrained channel [2]. In these cases, the 
required channel bandwidth is quantified by the number of bits being transmitted between the sensors 
and the FC. In contrast, the analog method consists of amplifying and then forwarding the sensed data to 
the FC, while respecting a power constraint O, JH. The transmissions can be appropriately pulse-shaped 
and amplitude modulated to consume finite bandwidth. The channels between the sensors and the FC can 
be orthogonal, in which case, the transmissions from each sensor are separately received at the FC f2\. 
On the other hand, with multiple-access channels between the sensors and the FC, the noisy sum of all 
the transmissions are received at the FC to make a decision (Si, O-IISl. The bandwidth requirements of 
sensor networks with orthogonal channels scale linearly with the number of sensors, whereas, when the 
channels are multiple-access, transmissions are simultaneous and in the same frequency band, keeping the 
utilized bandwidth independent of the number of sensors in the network. 

Distributed detection problems have been mainly studied assuming a single receive antenna at the FC. 
It is possible that introducing multiple antennas at a receiver may overcome the degradations caused by 
multi-path fading and noise. Inspired by conventional MIMO systems, a natural question is how much 
performance gain can be expected from adding multiple antennas at the FC in a distributed detection 
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spanias} @ asu.edu). 



2 



D 

e 
c 

s 



n 



Fig. 1. System Model: A random parameter is sensed by L sensors. Each sensor transmits amplified observations over fading multiple 
access channels to a fusion center with TV antennas. 

problem. However, this question cannot be directly answered by the studies in the MIMO literature. 
Adding multiple antennas to the FC for distributed detection problems is different when compared to the 
analysis of conventional MIMO systems for two reasons: (i) the presence of sensing noise (the parameter 
of interest is corrupted before transmission); and (ii) a large number of sensors enable asymptotic analysis. 
To the best of our knowledge, references flU, [[lO'l are the only works that consider multiple antennas at 
the FC. In dH, a decision fusion problem with binary symmetric channels between the users and the FC is 
considered where the data are quantized at the sensors, transmitted over parallel channels, and processed 
after being received by three antennas. In fTO], the authors consider multiple antennas at the FC. However, 
they consider a set of deterministic gains for the orthogonal channels, known at the sensors. They do not 
consider multiple-access channels, or characterize the performance benefits of adding antennas at the FC 
in the presence of fading. The system models in lfTn - lfT6ll are similar to adding multiple antennas at 
the FC, where the authors consider other forms of diversity, such as independent frequencies, CDMA 
codewords or several time intervals over fast-time- varying channels. The main difference between these 
papers and our results is the fact that we use asymptotic techniques to investigate the benefits of adding 
multiple antennas at the fusion center, when the number of sensors grows large. We show that the gain on 
the error exponent by adding antennas to the FC when there is no CSI at the sensors grows linearly with 
the number of antennas. In stark contrast, when there is CSI at the sensors, only limited gains are possible 
by adding antennas at the FC. This is unlike what we see in traditional MIMO wireless communications, 
where adding antennas at the FC will result either in diversity gain or array gain, for asymptotically large 
SNRs. 

In this paper, a distributed detection problem over a multiple access channel, where the FC has multiple 
antennas is considered (Figure [I])- The data collected by the sensors are transmitted to the FC using 
the amplify and forward scheme, with a total power constraint on the sensor gains. Performance is 
evaluated when the sensors have no channel information, have full channel information and partial channel 
information in the presence of fading, both with zero and non-zero mean. Analysis is performed for two 
cases: (a) large number of sensors and a fixed number of antennas, and (b) large number of antennas and 
sensors with a fixed ratio. In each case, the error exponent is used as the metric to quantify performance 
through the effect of channel statistics and the number of antennas. It is shown that the system performance 
depends on the channel distribution through its first and second order moments. This information is used 
to address our main objective, which is to quantify the gain possible by adding multiple antennas at the 
FC over fading multiple-access channels for distributed detection problems. 

II. System Model 

A sensor network, illustrated in Figure [H consisting of L sensors and a fusion center with antennas 
is considered. The sensors are used to observe a parameter G {0, 6}. The value, xu observed at the 
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sensor is 

under Hq 




xi = < (1) 
under Hi 

for / = 1,...,L. It is assumed that r]i ~ CA/'(0, cr^) are iid, the hypothesis Hi occurs with a priori 
probability, < pi < 1, and the hypothesis Hq with probability po = 1 — pi. The l^^ sensor applies a 
complex gain, ai, to the observed value, xi. This amplified signal is transmitted from sensor / to antenna 
n over a fading channel, hni, n = 1, A^, and / = 1, L, which are iid and satisfy -Efl/in^p] = 1. Unless 
otherwise specified, no other assumptions are made on the channel distribution. The n^^ antenna receives 
a superposition of all sensor transmissions in the presence of iid channel noise, z/„ ~ CJ\f{0,al), such 
that 

L 

yn = ^ hinai{Q + r]i) + (2) 

i=l 

where {r^j}^^ and {yn}n=i independent. 

Defining o; as an L x 1 vector containing {ajjfli, D(q;) an Lx L diagonal matrix with the components 
of a along the diagonal, the received signal is expressed in vector form as 

y = Hae + HD(a)r7 + u, (3) 

where H is an x L matrix containing the elements hni in the n^^ row and l^^ column, r/ is an L x 1 
vector containing {r]i}f'^i, and is an x 1 vector containing {h'n}n=i- Based on the received signal, y 
(from ([3])), the FC decides on one of the two hypotheses Hq or Hi. Since the FC has full knowledge of 
H and a, y is Gaussian distributed under both hypotheses: 

i/o:y~CA/'(0^,R(a)) 

Hi-.y^ CJ\f{eYl(x, R(q:)) (4) 

where Oa? is an x 1 vector of zeros and R(q:) is the N x N covariance matrix of the received signal 
given by 

R(a) = HD(q:)D(q:)^H^ + al I^. (5) 
We consider detection at a single snapshot in time, and therefore, we do not have a time index. 



A. Power Constraint 

The i^^ sensor transmits ai{Q + rji). The total transmitted power is given by 



Pt = E 



1=1 



a,- 



(6) 



1=1 



It should also be noted here that the instantaneous transmit power from the sensors is \ai{0 + rii)\'^ . This is 
a function of the actual realizations of sensing noise, making it difficult to predict and constrain. Therefore, 
we constrain a^'s, which allows imposing an average (over sensing noise) power constraint. The sensor 
gains, are constrained by 

Pt 



1=1 



Pi9^ + <yl' 



(7) 



4 



B. The Detection Algorithm and its Performance 

Given the received data, y, the FC selects the appropriate hypothesis according to 

^{dy^Ii{cx)-^Yi(x} % -d'^a.^Yi^Ii{a.)-^Yla. + r, (8) 

Ho 2 

where r is a threshold that can be selected using the Neyman-Pearson or the Bayesian approach. Using 
(HI) and ([8]), and the Bayesian test with the detection threshold, r = (1/2) ln(]9o/pi), the probability of 
error conditioned on the channel can be calculated as 

Pe\ii{N) = PoQ (to + r/co) + piQ {cu - r/cu) , (9) 



where uj ■.= 6y' a.^li^'R{a.) ^Hq:/2 for brevity, is the number of antennas at the FC and Q{x) = 
■fx" tI^^"^ ^"^dy. The error exponent is defined in terms of the conditional error probability for the FC 
with N antennas as [[Bl, IfTTll. ifTSl 

S{N)= lim -llogP,|H(iV). (10) 

Note that even though Pe\H.{N) in (|9l) is a channel-dependent random variable, we will show that the limit 
in (flOl) converges in probability to a deterministic constant for the cases of interest to us. Substituting ^ 
into (fTOl) , using L'Hopital's rule, and the Leibniz Integral rule for differentiating under the integral sign, 

S{N)= lim --e'^cx"U.^Ii{cx)-^licx (11) 

L— >oo 8 L 

in probability, which does not depend on pq and pi. Since S{N) is the negative exponent of the probability 
of error, a larger value represents better performance. The error exponent in (fTT)) is a deterministic 
performance metric over fading channels and depends on fading statistics. It can also be viewed as 
a "generalized SNR" expression in this system with multiple sensor and channel noise sources. We 
follow [[T5||. ifm . ifTSl in our definition of the error exponent in (flOl) . Alternatively, one can consider 
the unconditional error exponent, E-H{Pe,\-ii{N)], which would depend on the distribution of H in (flOl) . in 
place of Pe\ii{N). We will not pursue this approach herein. 

Our primary focus throughout this paper is the dependence of (fTTI) on 

(i) the number of antennas, N, for different fading-channel distributions; 

(ii) different assumptions about the dependence of the sensor gains, ol, on the channel, H. 

With the Neyman-Pearson test, rather than the Bayesian test, it can be shown that the error exponent 
is given by limi_j.oo 0.5L~^6'^Q;^H^R(a)~^HQ:, which does not depend on the false alarm probability 
and is a factor of four greater than the error exponent derived in the Bayesian case. Since the two cases 
differ only by a fixed constant, the Bayesian approach will be used throughout. 

III. Performance over AWGN channels 

The error exponent with AWGN channels is computed to establish a benchmark for the fading case of 
the next section, which is our main focus. For AWGN channels, = 1. Due to symmetry and to respect 
the power constraint, = -^/P/L, Vz. Defining 1^ as an L x 1 vector of ones, and Iatxl as an x L 
matrix of ones, we have a = \fPjL 1l and H = Ia^xl- Substituting these in ([5]), 



R:=R(/p7l1l) = cr'},PlNxN + ct^Iat. (12) 

The inverse of (fT2l) can be expressed using the Sherman-Morrison-Woodbury formula for matrix inversion 
and substituted into (fTTj) to yield 

gAWGN(iV):=^ ,, (13) 
8 A^7c + Pi7s + 1 
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where the sensing SNR is defined as 7s:=6'^/cr^, and the channel SNR, 'yc'-= Pt/ (^l- Since the partial 
derivative dSp^Q^{N) / dN > 0, for the AWGN case, having multiple antennas improves the error 
exponent which can be interpreted as array gain on the channel SNR 7c. As a special case, consider 
= 1, to get the result for the single antenna case: 

Wn(1) = I J""^" ^ , . (14) 
8 7c + Pils + 1 

With pi = 0.5, 7c = 1 and 7s = 1, adding a second antenna at the FC provides a gain of 3.1dB. Adding 
a third antenna provides a further gain of 1.34dB, indicating diminishing returns. To study the benefits 
of having multiple antennas, we compare the error exponent in each case with £^awgn(1)- The multiple 
antenna gain for the AWGN case is given by 

gAWGN(iV) N^c + iVpi7s + N 
<-^AWGNUVj := — — = — ■ — — . (n) 

cawgn ( 1 ) A^7c + Pi7s + 1 

It can be seen from (fT5l) that by making sufficiently large, and 7c sufficiently small, (fT5l) can be made 
arbitrarily large. In contrast, it will be seen in Section IIV-BI that when the channels are fading and known 
at the sensors, the corresponding gain expression will be bounded for all parameter values, indicating 
limited gains due to antennas. 

IV. Performance over Fading Channels 



Suppose that the elements of H are non-zero-mean: = ^jKjiK + 1) + + \)h^f , where 

the first term is the line-of-sight (LOS) component, h^f is the zero-mean diffuse component, and the 
parameter K is the ratio of the LOS power to the power of the diffuse component, chosen so that the 
channel satisfies E^Ki(^\ = S[Kf P] = 1. 

In what follows, different cases of channel state information at the sensors (CSIS) are considered. 

A. No Channel State Information at the Sensors 

When the sensors have no channel knowledge, then the sensor gains are set to a = PjLXi^ due to 
the i.i.d. nature of the channels and to respect the power constraint in Substituting in ([5]), 



R := R(/p7l1l) = olp\w^ + a^I^. 



(16) 



L 

Since the elements of H are i.i.d., from the weak law of large numbers, 

lim R = alP——lN^N + ^ '-In, (17) 

L-5-00 ' A + 1 A + 1 

in probability. Since the right-hand-side of (flTl) is non-singular, limL_>oo R~^ = (liniL^.oo R)~^ [fT9l Thm. 
2.3.4]. Using the matrix inversion lemma on (flTI) and substituting into (fTTI) . 

P (N K^ PjK+l) , >f l>f - 



1=1 



9^ alP^K{K + 1) 

lim 

8 [a2p + al{K + 1)] [a^^PNK + a^P + al{K + 1)] 



^ N L 
n=l 1=1 



(18) 



Using the weak law of large numbers and ([7]), the error exponent can be expressed in terms of 7c and 7s 
as 

8 7c(A^A: + 1) + (pi7s + 1) (AT + 1) 
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which can be shown to be a monotonically increasing function of A^, K, % and 7c, as expected. For 
the single antenna case, using (fT4l) we have £^noCSIs(1; = ^^awgn(1)-^/(-^ + 1): which is a factor 
K/{K + 1) worse than £^awgn(1)- 

As the antennas increase, liniiv^oo ^^NoCSIs(^, = 7s/8, which is the same as limAr_j.oo i^AWGN(A^)- 
That is, so long as there is some non-zero LOS component, as the number of antennas at the FC 
increases, the performance approaches the AWGN performance even in the absence of CSI at the sensors. 
Furthermore, it can be seen that lim;^_j.oo ^NoCSIs(^) K) = £awgn{N), which matches the AWGN result, 
as expected. 

To characterize the gain due to having multiple antennas at the FC, we define 

r '-r T.^ _ £noCSIs{N,K) iV(K + l)(7e + pi7s + l) .... 

L^NoCSis U V , A j . - — - — rTuzFVT\ ■ '^'^^^ 

^^NoCSis(l, K) -ic[NK + 1) + (pi7s + 1){K + 1) 

When the channel noise is large, (7c — ?■ 0), we have GnoCSIsI^; K) = N and the gain increases with the 
number of antennas at the FC. However, when 7c — t- 0, the absolute performance of the system is poor, as 
can be verified by substituting in (fT9l) . Conversely, when the channel SNR grows, the maximum gain in 
(|20|) is given by (/^ + I)/ K. This leads to the conclusion that when the channels between the sensors and 
the FC are relatively noise-free, there is little advantage in having multiple antennas at the FC when K is 
large. When the channel is zero-mean {K = 0), the error exponent in (fT9l ) is zero for any A^, indicating 
that the probability of error does not decrease exponentially with L for any A^, confirming results from 
S, [|20l , [|2n . However, from (|20l ). it is clear that the gain satisfies lim^^-.^o G'noCSIsI-^! = which 
shows that when the channel is zero-mean, gain in the error exponent due to antennas is linear and can 
be made arbitrarily large. We have thus established the following: 

Theorem 1: For zero-mean channels, with no CSI at the sensors, the error exponent in (fT9l) is zero and 
therefore, the error probability does not decrease exponentially with L for any number of antennas, A^. 
The antenna gain, defined in (|20l ) satisfies limx^o G'noCSIs(^) = N, implying unlimited gains from 
multiple antennas for zero-mean channels when CSI is unavailable at the sensors. 

In what follows, it will be seen that when CSI is available at the sensors, the antenna gain is bounded 
over all parameter values for zero-mean channels. 



B. Channel State Information at the Sensors 

We have just seen that when the non-zero-mean channel assumption does not hold, the incoherent 
sum of signals at each each antenna leads to poor performance at the FC, which results in a zero error 
exponent. If channel information is available at the sensors, the sensor gains can be adjusted in such a 
way that the signals are combined coherently. It should be noted here that full CSI at the sensors implies 
full CSI of the network, H, at the sensors. In such a case, ol is chosen as a function of the channels, H. 

As a benchmark result for fading channels, the sensor gains are selected in such a way as to maximize 
the error exponent of the system given in (fTTI) . subject to the power constraint in ([7]): 

ckqpt = argmax [(x^H^R{a)-^Ha] subject to ||af < P, (21) 
to obtain the error exponent in the presence of CSIS, 

^^csis(iV) = lim ^yttgpTH^ (a2HD(aoPT)D(aoPT)''H^ + ^^Ijv)"' HaoPT- (22) 

The optimization problem in (|2TI) is not tractable when A^ > 1 since R(q;) depends on H and a. In 
order to assess the effect of number of antennas, the solution for (l22b with A^ = 1, and two upper bounds 
on (|22|) are derived for A^ > 1. 
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1 ) Solution for Single Antenna at the FC: When = 1, the channel matrix reduces to a column vector, 
given by [hih2 . . . h^]^, where hi is the channel between the i-th sensor and the FC. The maximization 
problem in (|2TI) reduces to 

2 



ckqpt = argmax ■ 



i=l 



i=l 



subject to 



a,; 



< p. 



(23) 



\hi\ +cr„ 



i=l 



A similar problem was formulated in [|7]| and in a distributed estimation framework in [|3]|, [[22l. We 
recognize that the best value for the phase of the sensor gain is Zai = —ipi where tpi = Zhi. Therefore, 
we set Zai = —ipi,\/l. We then define s := J2i=i ^i^i swap the objective function with the constraint 
so we can rewrite the optimization problem as 

L L 



olqpt = argmin^^ 



subject to cr^ \ai\^ \hi\^ + 1 < VtS^ 
1=1 

L 



0, 



(24) 



1=1 



where Vt is an auxiliary variable. The optimization problem in (|24l) is now a (convex) second-order-cone 
problem ll23l . Using the Karush-Kuhn-Tucker conditions ll23l . the optimal solution is given by 



P 



E 



I ft, 





\K 









(25) 



^Z^^\p\h,\^a^^ + al 
The error exponent can be obtained by substituting (l25l) in (|22)) with = 1: 



8. 



CSIS 



lim ^ 5 — 

1=1 ^r, ^ p|/i,|2 



E 



(26) 



from the weak law of large numbers, where the expectation is with respect to {hi}. As an example, for 
Rayleigh fading channels ^ yields JM §3.353] 



^^csis(l) 



1 

32 



7s 



' Plls + 1 

2 exp 

7c 



Pi7s + 1 
27c 



El 



Plls + 1 

27c 



(27) 



where -Ei(-) is an exponential integral function [25", pp. 228]. The expression for £^csis(l) is obtained 
when the channels between the sensors and the FC are fading. To compare with the AWGN case, note that 
Px/ {(j'^Px + al) in (|26l) is a concave function of x, and from Jensen's inequality, £^awgn(1) > ^^csis(l), 
as expected. 

Since (|26l ) is rather complicated, it is desirable to find a simpler expression as a lower bound to (|26l ). 
Any choice of ||q;||^ = P will yield such a lower bound, since ckqpt is optimal. Considering phase-only 
correction at the sensors, at = \J P/Lexp{—jZhi) is substituted in (fTT)) with = 1 to yield the error 
exponent for phase-only CSIS for = 1: 

2 



P 



s 



PO\ 



lim 



9' 



cr; 



(28) 



1=1 
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From the weak law of large numbers, the random sequences in the numerator and denominator converge 
separately. However, since the expression for £po{l) is a continuous function of these sequences, the 
value of <£^po(l) converges to [|26i, Thm. C.l] 

Spoil) = {E[\hi\]fSAWGM (29) 
in probability, since -Efl/i/p] = 1. The expression in (|29l) serves as a lower bound to £^csis(l) as follows: 

^^awgn(1) < ^csis(l) < ^awgn(I), (30) 

where C = (^[|/i/|])~'. 

2) Upper Bound (AWGN channels): Since (|2T]) cannot be solved in closed form when > 1, one 
cannot evaluate the error exponent in (|22)) by substitution as it was done for = 1. Two upper bounds on 
((22)) will be convenient at this stage. Since the AWGN performance is a benchmark for fading channels, 
the error exponent of the system over AWGN channels is an upper bound on that of fading channels, 
even in the case of full CSIS. Therefore, the first upper bound to (|22l) is given in (fT3l) : 

^csis(iV) < WN(iV) = \ , , ^J"^" ^ , . (31) 



3) Upper Bound (No Sensing Noise): Clearly, (|22)) is a monotonically decreasing function of the sensing 
r^. The second benchmark is obtained by setting cr^ 



noise variance, a^. The second benchmark is obtained by setting cr^ = 0, which also affects aopT in (|2T)) . 



since R(q;) no longer depends on a when cr^ = 0. Substituting this in ((2T)) . the optimal value of a when 

argmax (a^H^Ha) subject to ||q:||^ < P. (32) 



= is 



The solution to (|32l) is the eigenvector corresponding to the maximum eigenvalue of H^H, scaled in a 
way to satisfy the constraint with equality. Substituting into (l22l) with cr^ = 0, we have the second upper 
bound to i^csis(^): 

5(iV, K) = ^-^ lim A^ax fTH^'H^ , (33) 



----- - 'iiiax 1 T 

8 af, L^oo \L 

where Amax(-) denotes the maximum eigenvalue function. Using the fact that Amax(H^H) = Amax(HH^), 
and that Ainax(-) is a continuous function of the matrix elements [fT9i Thm. 8.1.5], one can interchange 
the limit with the maximum eigenvalue function [26, pp. 422, Thm. C.l] to yield 

B{N, K) = ^4^-- ( lim yHH^) . (34) 

From the weak law of large numbers, 

1 K 1 

lim — HH = — — -IjvxAT + ^T— riAfxTV, (35) 
L-s>oo L K + 1 K + 1 

in probability, so that with the substitutions cr^ = and O'^P/a'l = 7c/pi, we have the bound: 

1 ^ NK + 1 

^csis(iV) < B{N, K) = --^r^. (36) 

8pi K + 1 

In (|36l) , B{N, K) is an upper bound when there is sensing noise in the system. When there is no sensing 
noise, it is the actual error exponent of the system with full CSIS. Furthermore, Yim.K^ooB{N,K) = 
lim^^_>oo ^^awgnI^)^ verifying that as K ^ oo, B(N,K) converges to the AWGN error exponent with 
no sensing noise. In addition, if i^' = 0, there is no advantage to having multiple antennas at the FC, for 
asymptotically large number of sensors, since the right hand side of (l36l) is independent of in that case. 
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Since both Sawq-^{N) and B{N, K) are upper bounds to i^csis(^): ^ combination of the two bounds 
T[vm.[8p^Q^{N)^B{N,K)], provides a single, tighter upper bound. Equating the right hand sides of (1311 
and (|36l) . it can be shown that this combined upper bound is given by 



B[^,K} If (^r,<N{NKTr) 



2 \ N-1 



Combining the upper and lower bounds, 

1 



^SawgM < ^csis(l) < ^csis(iV) < C(iV, K), (38) 



obtained from (|27I) . (1301) and (|37I) . The bounds in (1381) will be used to further examine the effect of on 

^csis(iV). 

The value of <£^csis(^) from (|22|) is the best achievable performance for fading channels. Defining the 
gain due to multiple antennas as G'csis(^) := ^csis(^)/^csis(l)5 the following theorem can be stated: 

Theorem 2: When the channels have full CSI at the sensors, the gain due to multiple antennas at the 
FC can be upper bounded as 



Gcsis(A^) < Ct. 7TT < Cmm 

t^AWGNl-Lj 



Nz + l'^^ ' K + l 



(39) 



where 2;:=7c/(pi7s + !)• 

Proof: The first inequality in (|39| ) follows from the first inequality in (l38l) . The second inequality in 
(|39| ) follows from the last inequality in (l38l) and dividing the terms of (l37l) by (fT4l) . ■ 

With j9i = 0.5, K = 1, 7^ = 1 and 7s = 1, for = 2, ^0313(2) < 1.4286C. For iV = 3, Gcsis(3) < 
1.6667C and for = 4, Gcsis(4) < 1.8182C. These results indicate that there is diminishing returns in 
the multiple antenna gain. 

Corollary 1: G'csis(iV) can be bounded by an expression depending on and K only: 



Proof: The first argument of the min[-, ■] function of the right hand side of (1391) is a decreasing 
function in z and the second argument is an increasing function in z. Therefore, when the arguments are 
equal for fixed values of A^ and K, the maximum value of the min[-, •] function is obtained. This occurs 
when z = N^^{NK + 1)~^(A^ — 1), allowing us to upper bound the min[-, ■] function by the value in 
(HOI). ■ 
Corollary 2: When the channels have zero-mean, the maximum gain due to having multiple antennas 
at the FC is bounded by a constant independent of A^ and only dependent on ^ = (i?[|/ii|])^^: 

GcsisiN) < 2C. (41) 



Proof: Substituting i^' = 0, it is clear that (|40|) is monotonically increasing in A^. Taking the limit as 
A^ — 7- 00 yields the proof. ■ 

As an example, in the case of Rayleigh fading, when full channel information is available at the sensors, 
the maximum gain that can be obtained by adding any number of antennas at the FC for any channel or 
sensing SNR is at most 2( = S/ir, which is less than 3. 

The results in (|39l)- (|4T1) have been derived for the case of iid sensing noise. We now address the 
correlated sensing noise case. To this end, we define as the L x L covariance matrix of the sensing 
noise samples, {r]i}[^^. 
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Theorem 3: Suppose that the sensing noise samples are correlated and let Amin be the minimum 
eigenvalue of R^. The gain due to multiple antennas in (|39l) holds with the change 5; = 7c/(pi7s + 1), 
where 7s := 6'VAmin- 

Proof: The proof is shown in Appendix |Al ■ 
Theorem [3] shows that any full-rank sensing noise covariance matrix changes the conclusion in (|39l ) 
only through a redefinition of z. By maximizing over z, the same upper-bound in (l40l) is obtained, and 
for zero-mean channels, the bound in (|4TI) remains valid. This shows that the bounds in (l40l) and (|4TI) 
are general, and hold even when the iid condition is relaxed to any arbitrary full-rank covariance matrix, 
R^. The gain due to adding multiple antennas is still upper-bounded by a factor of 2C,, for zero-mean 
channels, when there is full CSI at the sensors. 



C. Phase-only CSIS 

One simplification to the full CSIS case is to provide only channel phase information to the sensors. 
For the single antenna case, and when the channels between the sensors and the FC have zero-mean, the 
phase-only results have been presented in (|29l ) and (|30l) . What follows is an extension of those results to 
the multiple antenna case when K = 0. 

Since there is only phase information at the sensors, the amplitudes of the sensor gains are selected 
such that \ai\ = ^P/L, V/, so that D(a)D(cK)^ = {P/L)Il and R( a) is given by (fT6l) . 

With phase-only information, one can constrain \ai\ to be constant to reformulate (|2TI) as the following: 

p 

CKpo = argmax a^H'^^Ha subject to |aj|^ = — , z = 1, 2, . . . , L. (42) 
In Section |IV-E4[ a semidefinite relaxation approach will be presented to solve (|42l) . 



D. Asymptotically large sensors and antennas 

When CSIS is available, (l39l - |4T]) shows that only limited multiple antenna gains are available. It 
is interesting to see whether such limits would still be present if — )■ 00 simultaneously with L. A 
similar problem was considered, but in the context of CDMA transmissions in [1I6I . Note that this will 
in general yield results different than first sending L — t- 00 and then — )■ 00 as was done in Section 
IIV-BI Such a situation can be interpreted as a case where a group of sensors is transmitting to another 
group, functioning as a virtual antenna array ||27l . For such a system the scaling laws when L and A^ 
simultaneously increase [|28l pp. 7], in such a way that 

lim 4 = /3, (43) 

are of interest. It should be noted that in spite of scaling the number of sensors and antennas, the power 
constraint is still maintained. 

In this case, the error exponent is redefined as 

£°°(/3)= lim -llogP,|H(Ar), (44) 

L,Af— 5>oo L 



with (1431) satisfied. Similar to the upper bounds in (|3TI) and (1361) . upper bounds on (|44l) are now derived. 
For the AWGN case, 

^°°(/3)<^:rwGN:= lim WnI^V) = lim \ ^^'^^ = ^7s. (45) 

L,Af-s>oo L,Ar-s>oo 8 A7c + Pi7s + 1 8 

When there is no sensing noise, with af^ = 0, the second bound can be calculated as 

S^{P)<B^{P):= \im ^:^A,„ax f jH^hV (46) 
L,N^oo af, \ L J 
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For fading channels with i^' > 0, it can be shown that the error exponent in (|46l) goes to infinity. Therefore, 
with any line-of-sight (LOS) and no sensing noise, increasing the number of sensors and the number of 
antennas to infinity provides very good performance. When K = 0, the Marcenko-Pastur Law [i28, pp. 
56] provides an empirical distribution of the eigenvalues of A^^H^H. From [i29l , [|30l . the maximum 
eigenvalue of A^^^H^H is shown to converge in such a way that 



lim Amax 



H 



N \VN 



H 



2 



(47) 



in probability, which yields 



= 12£<1±/^. (48) 

8 Pi l3 

which is the optimum performance of the system in the absence of sensing noise. Similar to (|37l) . the 
minimum of (l45l) and (l48l) yields 

bs if Pa,^>^^ 



8pi /3 -^^'J - (l+v^)2 



The gain due to antennas is expressed in terms of the ratio /3 in (l43l) as G°^(/3) : = i^°^(/3)/£csis(l)- Using 
the bounds, we have the following: 

Theorem 4: With asymptotically large number of sensors and antennas, the gain due to having multiple 
antennas at the FC is bounded by 

G°°(/3)<c(^l+ ^^^/^^ (50) 

Proof: The relationship between £^awgn(1) and i^csis(l) from (|30l ) provides a lower bound on 
^0313(1)5 and consequently an upper bound on G^{f3), to yield the first inequality in (|5T1) below. The 
expression in (|49l) provides an upper bound on £°^(/3), and dividing by (fT4l) yields the second inequality 
in 



£°°f/?) 

G°^(/3)<C p <Cmin 

t^AWGNUJ 



(51) 



where w:=7c/(pi7s + 1). The first argument in the min[-, ■] function decreases as w increases, while 
the second argument is an increasing function of w. Therefore, the min[-, ■] function is maximized when 
arguments of the min[-, ■] function are equal for a fixed value of /5. This result is obtained when w = 
(1 + v^)-2/3, to yield ^ and the proof. ■ 
To interpret (l50l) . cases corresponding to three values of (3, are considered: 

(i) /9 -C 1 (A^ scales faster than L): When the number of antennas increases at a faster rate than 
the number of sensors, it can be seen that B°°{(3) is large. When there is no sensing noise, the 
performance obtained is exactly B°°{/3) as seen in (l48l) . In this case, arbitrarily large gains are 
achievable. In case there is sensing noise in the system, ^^^qn and B°°((3) become bounds, and the 
gain is bounded as shown in (l50l) . As /3 — )■ in this case, the bound goes to infinity, which indicates 
that there could be large gains possible. 

(ii) (3 = 1 {N scales as fast as L): The number of antennas at the FC and the number of sensors scale 
at the same rate, the maximum possible gain can be calculated from (|50l) to yield G°^{1) < 5^- 

(iii) /3 1 (A^ scales slower than L): When the number of sensors scales much faster than the number 
of antennas at the FC, it resembles the previous setting where L — )• 00, first, and A^ was scaled. Not 
surprisingly, when /3 is large in this case, G°^{I3) < 2(, same as in Section ITV-BI 





GnoCSis{N,K) from m 


Gcsis(N,K) from m 


G°°(/3) from 150) 


K>Q 


0{N) when 7c = 0; 0(1) when 7c > 


0{N) 


Undefined 


K ^0 


0{N) 


0(1) 


0(/3~") as /3 ^ 0; 0(1) as /3 ^ oo 



TABLE I 

Order of gain due to multiple antennas at the FC for large number of sensors, L. 



It should be noted here that in cases (ii) and (iii), where both the number of sensors and antennas are 
scaled to infinity simultaneously, only limited gain is achievable, when the sensors have complete channel 
knowledge. 

In Table |I] we summarize the rate at which the gain due to number of antennas increases, both when 
CSI is available and unavailable at the sensor side. Recalling that the gain is defined in terms of the ratio 
of error exponents relative to the single antenna case, all the results in the table apply when L is large, 
which is a major distinguishing factor between this study and standard analysis of multi-antenna systems. 
It is seen that when K > the gain in error exponent grows like 0{N) depending on whether CSIS 
is available and whether 7c = 0. More interestingly, when the channel is zero-mean (K — )■ 0), adding 
antennas improves the error exponent linearly when CSIS is not available. In stark contrast, when CSIS 
is available, the gain is bounded (0(1)) by 2(^. Finally, the column on the right of Table IJ illustrates how 
the gain depends on the ratio (3 = L/N as both and L increase. The error exponents for K > are 
infinite, yielding an undefined gain. For zero-mean channels, the dependence on /3 indicates an increasing 
gain when /3 is small (L <^ N), and bounded gain when /3 is large (L ^ N). 



E. Realizable Schemes 

So far, we have provided bounds on the achievable gains due to antennas when CSI is available at the 
sensors, without providing a realizable scheme. This is because the calculation of aopx in (|2TI) in closed 
form is intractable. Moreover, it is not clear how ex. should be chosen as a function of H when > 1 to 
achieve a multiple-antenna gain. This is because each sensor sees channel coefficients, corresponding 
to antennas, and each channel coefficient has a different phase making the choices of Za^ non-trivial. 
We now present two sub-optimal schemes for the full CSIS case that are shown to provide gains over the 
single antenna case. 

1) Method I: Optimizing Gains to Match the Best Antenna: In this method, the sensor gains, ex., are 
selected in order to target the best receive antenna. However, the received signals at all of the other 
antennas are also combined at the FC, which uses the detection rule defined in ([8]). 

Since L is finite for any practical scheme, (l25l) will be used to select ex and (|26l) without the limit can 
be used to assess which antenna has the "best" channel coefficients. Therefore, using the channels from 
the sensors to all of the receive antennas, 

n* = argmax 5 — , (52) 

1=1 + PK7F 

is calculated and the sensor gains are set to (l25l) computed for the channels {hnH}i=i- The FC then uses 
all of the receive antennas for detection using ([8]). Since there are multiple antennas at the FC, for any 
realization of the channels between the sensors and the FC, the error exponent of this scheme is at least 
as good as the single antenna case. 

Such an approach requires the calculation of (|52|) and the corresponding at from (|25l) . Since these 
calculations require the complete knowledge of H, they can be calculated at the FC, and fed back to the 
sensors. 
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2) Method II: Maximum Singular Value of the Channel Matrix: It was shown in Section IIV-B3I that 
when (7^ = 0, the bound obtained in (|36l ) is achievable. In this method, the values of ex. are selected 
as though there is no sensing noise. The sensor gains, a, are selected in such a way that they are a 
scaled version of the eigenvector corresponding to Amax (H^H), such that ||q:P = P. In most practical 
cases, sensing noise is non-zero, and therefore, this method is sub-optimal. Similar to Method I, ex. can 
be calculated at the FC and fed back to the sensors. 

3) Hybrid of Methods I and II: Since Method II is tuned to perform optimally when there is no sensing 
noise, it outperforms Method I when the sensing SNR, 7s, is high. As the sensing SNR reduces. Method 
I begins to outperform Method II. These observations are illustrated and elaborated on in the simulations 
section (Section |Vl Figure [8]). 

Since one of the schemes performs better than the other based on the value of 7s, a hybrid scheme 
can be used: Method I for low values of 7s, and Method II for high values. The exact value where the 
cross-over occurs depends on the parameters of the system, and can determined empirically. An example 
is shown in the simulation section in Figure [8l where it is also argued that an underestimation of the 
value of 7s is tolerable, while an overestimation is not. 

4) Semidefinite Relaxation: Following |[22l . [|3T1l a semidefinite relaxation of the problem in (l42l) is 
obtained as follows: 

Xpo = argmaxtrace(H'^HX) subject to X ^ 0, 

X 

_P , _ 

Xjj = — , i = 1, 2, . . . , L, (53) 

where X is an L x L matrix. If X has a rank-1 decomposition, X := olol^ , then a is a solution to (|42|) 
[|22]| . [|3T1 . In the more likely case where X does not have rank-1, then an approximation to the solution 
of (|42)) is obtained by choosing cx as the vector consisting of the phases of the eigenvector corresponding 
to the maximum eigenvalue of X. The semidefinite relaxation in (l53l) causes a loss of upto a factor of 
7r/4 in the final answer of (|42l) [[311 . The phases of eigenvector corresponding to the maximum eigenvalue 
of Xpo are extracted to constitute a possible set of values of a. In order to obtain the solution to the 
SDR problem, an eigenvalue decomposition of Xqpt is required, which is an 0{L^) operation [19J. It is 
argued with the help of simulations (Figure (9]) that the SDR outperforms the hybrid scheme when 7s is 
small, at the expense of increased complexity. 

V. Simulation Results 

The theoretical results obtained are verified using simulations. The channels are generated as complex 
Gaussian (Rayleigh or Ricean) for the purposes of simulation, even though the results only depend on 
the first and second order moments of the channels. 

In Figure [21 it is verified that increasing the number of sensors improves the performance except when 
the channels are Rayleigh fading and there is no CSIS. Since the error exponent is zero for the Rayleigh 
fading case with no CSIS, the asymptotic average probability of error is computed and plotted. The Ricean 
case outperforms the Rayleigh fading case, and the AWGN channels provide the best performance. It can 
also be seen that the decay in probability of error is exponential in L, when the channels between the 
sensors and the FC are AWGN or Ricean fading. The decay is slower than exponential when the channels 
are Rayleigh fading. This confirms the observations in Section [rV-A[ In all cases, the performance improves 
as the number of antennas increases. 

In Figure [3l the expression of error exponent is compared against the value of L^^ logPe|H(5) for 
increasing L, with AWGN channels and Ricean fading channels between the sensors and the FC. It can 
be seen that fewer than 200 sensors are required for the asymptotic results to hold. Therefore, in subsequent 
simulations, L = 200 sensors have been used. 

The effect of increasing the number of antennas on the error exponent for the AWGN case and Ricean 
fading case with no CSIS is seen in Figure SI As expected, increasing 7c improves performance and there 
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Fig. 2. Monte-Carlo Simulation: E[P^^^^yi(N)] for AWGN channels, Rayleigh fading channels and Ricean channels with no CSIS. 
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Fig. 3. Monte-Carlo simulation - Error exponent for AWGN and Ricean Fading channels. 



is an improvement in perfomiance as the number of antennas at the FC increases. As predicted in Section 
IIV-AI with an increase in N, the performance of £awgn{N) and i^NoCSis(^) get closer to each other. 
There is a large performance gain between the = 1 case and the N = 2 case, and almost the same gain 
between the N = 2 case and the A^ = 10 case, indicating diminishing returns, corroborating the results 
in Section Unl 

In Figure [51 the error exponent is evaluated when there is a single antenna at the FC. The cases of 
AWGN channels, Ricean channels with no CSIS, Rayleigh fading channels with full CSIS and Rayleigh 
fading channels with phase-only CSIS are compared in Figure [51 It is seen that the AWGN performance is 
the best, and when the Ricean channels have larger line of sight, the performance improves, as expected. 
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Fig. 4. Error exponent vs 7s for N = 1,2, 10 for AWGN channels and Ricean channels and no CSIS. 
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Fig. 5. Optimal Rayleigh performance, AWGN performance and Ricean no CSIS performance with one antenna at the FC. 



In fact, by increasing the amount of LOS, the no-CSIS Ricean case performs better than the full CSIS 
Rayleigh channel case, when 7c is large. The performance of the Ricean no CSIS case is a constant 
factor K/{K + 1) worse than the AWGN case, corroborating the result of f^NoCSisll, -^)- Similarly, the 
performance of the phase-only CSIS case confirms the result in (l30l) . For Rayleigh fading channels, the 
phase-only CSIS case performs a constant 7r/4 worse than the AWGN case. 

For the case of full CSIS, but with multiple antennas at the FC, bounds were derived on the error 
exponent of the system in Section IIV-B2I and Section IIV-B31 and combined to provide a single bound in 
(l37l) . The value of ^^csis(l) is set as a lower-bound on ^csis(-^)- In Figure [6l with N = 1, the upper 
bound can be seen to be about 0.76 dB (in terms of error exponent) away from the actual value at 7c = 8 
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Fig. 6. For a single antenna, optimal performance and performance bounds. 




Fig. 7. Comparison of antenna gains vs A*'. 



dB. For small values of 7c, the AWGN bound is better, and as 7c increases, the bound with the no sensing 
noise assumption is better, as expected. 

Figure |7] shows the effect of increasing the number of antennas at the FC on the antenna gains of the 
different systems. Also, for the cases of partial CSIS and full CSIS, the upper bounds on the antenna 
gains are plotted. The actual error exponent for the AWGN case is larger than for the Ricean no-CSIS 
case. However, as seen in Figure Ul the gain for the Ricean no-CSIS case is larger than the gain for the 
AWGN channel case. The bound on the Ricean CSIS antenna gain grows rapidly with N, as predicted by 
(|40l ). The maximum gains possible for the Rayleigh CSIS case and the Rayleigh no CSIS cases are also 
plotted. These results indicate that with full CSIS, there is not much to be gained by adding antennas at 
the FC, corroborating our results in Section IIV-BI 

The schemes introduced in Section IIV-EI for the known CSIS case are simulated in Figure [8l The 
performance of these schemes are evaluated for = 5 and = 50. The performance of these systems 
is compared against a lower bound given by £^csis(l) from (ITTI ) and an upper-bound, C(5, K) from (|37] ). 
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Fig. 8. Practical Schemes for iV = 5 and iV = 50 vs. fcsisCl) and C(5, 1). 

The hybrid scheme from Section IIV-E3I selects the better of the two practical methods depending on the 
value of 7s. It can be seen that even with these simple sub-optimal practical schemes, the hybrid scheme 
is always better than f^csisll)^ indicating that it is possible to obtain multiple antenna gain. However, 
for each A^, the hybrid scheme does not approach the upper-bound of C(5, K). When = 5, this is an 
expected result, since firstly, C (N, K) is a bound that is not necessarily achievable, and secondly, the 
practical schemes are obtained as sub-optimal approximations to the optimal scheme with full CSIS. The 
hybrid scheme for iV = 50 provides more gain over i^csis(l) than the hybrid scheme for = 5, but does 
not beat C{h,K). This means that although gains are possible with the practical schemes, large gains 
are not possible, as predicted by the bounds in Section IIV-BI For the hybrid scheme. Method I is better 
at low values of 7s and Method II is better at high values of 7s. The value of 7s at which the hybrid 
scheme changes methods can also be seen in the simulations. In Figure [8l the system has a channel SNR, 
7c = 10 dB, pi = 0.5 and the Ricean-fC parameter is one. When there are five antennas at the FC, 
the hybrid scheme changes from Method I to Method II at 7s ~ 3 dB, and when = 50, the change 
occurs at 7s ~ 8.25 dB. It can be seen that the hybrid scheme changes from Method I to Method II at 
different values of 7s based on the system parameters. It can also be seen that when Method I is selected 
by the hybrid scheme, the error in performance between Method I and Method II is small. However, 
when Method II is selected by the hybrid scheme, the performance gap between Method I and Method 
II increases rapidly as 7s increases. Therefore, an underestimation of the value of 7s is tolerable, while 
an overestimation is not. 

The semidefinite relaxation (SDR) approach in Section IIV-E4I is compared against the hybrid scheme 
(Section IIV-E3I ) in Fig [9l For the SDR solution, the value of Xqpt from (|53] ) is calculated using CVX, 
a package for specifying and solving convex programs in MATLAB [32|. It can be seen from these 
simulations that for low values of sensing SNR, 7s, the SDR solution outperforms the hybrid scheme. 
However, as the value of 7s begins to increase, the hybrid scheme (which is designed to be optimal 
as 7s — oo) outperforms the SDR solution. The comparison with the upper-bound on the optimal error 
exponent, C(A^, K) is tight with respect to the better of the hybrid and SDR approaches. In order to obtain 
the solution to the SDR problem, an eigenvalue decomposition of Xqpt is required, which is an 0{L^) 
operation [fT9ll . The SDR outperforms the hybrid scheme when 7s is small, at the expense of increased 
complexity. 
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Fig. 9. Hybrid realizable scheme, SDR relaxation and C{N, K) vs 7s. 



VI. Conclusions 

A distributed detection system with sensors transmitting observations to a fusion center with multiple 
antennas is considered. The error exponent is derived from the conditional probability of error. AWGN 
and fading channels are considered between the sensors and the fusion center with varying amounts of 
CSI at the sensors. A large number of sensors and a finite number of antennas, or a large number of 
sensors and antennas are considered. The gain due to having multiple antennas, defined in terms of the 
ratio of error exponents relative to the single antenna case for large L is quantified. The asymptotic nature 
of the sensors and the corruption of the sensed data with noise are the two major distinguishing factors 
in this work, compared to the conventional analysis of point-to-point MIMO systems. 

The gains due to the number of antennas at the FC, both when CSI is available and unavailable at 
the sensor side, are now summarized. It is seen that when K > when CSIS is unavailable, and when 
7c = 0, the gain in error exponent can grow unbounded with A^. When the channel is zero-mean (such 
as the Rayleigh case with K — 0), and when CSIS is not available, adding antennas improves the error 
exponent linearly. However, when CSIS is available, the gain is bounded by 2{E\\hi\]Y'^. In the special 
case of Rayleigh fading channels, when CSIS is available, the gain due to multiple antennas at the FC 
is bounded by S/vr. As both N and L increase, the error exponents for K > are infinite, yielding 
an undefined gain. For zero-mean channels the dependence on (3 from (|43l) indicates an increasing gain 
when (3 is small (L ^ A^), and bounded gain when (3 is large (L 3> A^). When the sensors use only 
channel phase information, the gain due to multiple antennas at the FC is upper bounded by (E'd/i;!])^^, 
irrespective of the number of antennas at the FC. 

In the case when full CSIS is available, it has been argued that having multiple antennas at the FC 
provides limited gains, which can be exploited by sub-optimal schemes. In one approach, the system is 
configured to beamform to the antenna that provides the best performance, where the FC still uses the 
data gathered at the other antennas. On an average, this is shown to perform better than in the single 
antenna case. Another approach is to assume there is no sensing noise, and set sensor gains tuned for 
such a system even when sensing noise is present. In this case, the system performs optimally when the 
sensing noise in the system is low. A hybrid scheme is proposed which selects the better of these two 
methods depending on the sensing SNR, 7s. 

Depending on the number of sensors and antennas at the FC, and their rates of growth, the following 
system design recommendation can be made. If CSIS is available and N <^ L, then for better performance, 
it is recommended to increase the number of sensors, rather than the number of antennas at the FC. 
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However, if can be increased at a much faster rate than L, it is possible to achieve greater gains due 
to adding antennas at the FC. 

Appendix A 
Proof of Theorem [3] 

We begin by noting that the presence of correlation in 77; affects the total average transmit power. 
Therefore, to prove Theorem [3l we need to reconsider the following in presence of correlation: (i) the 
power constraint; (ii) the AWGN upper-bound in (|3T1) ; (iii) the "no sensing noise" upper-bound in (|36l) . 
which will then be used to redefine the combined upper-bound in (|37] ). 

(i) Power constraint: The total transmitted power is given by 



E 



1=1 



{piOHl + R„) a, (54) 



and constrained as 
If ([55]) holds, then 



{pie^L + R^) a < Pt. (55) 

~ PlO^ + Amin ' ^^^^ 



also holds. Since (1561) is less stringent than (1551) . if (1561) is used instead of the original power constraint 
in (1551 ). an upper-bound will be obtained in the subsequent derivation of the error exponent. 

(ii) Upper-bound (AWGN channels): Recall that in this case, H = Iatxl- Since the sensing noise is 
not iid, a has to be selected in such a way that the error exponent is maximized: 

maximize 0!.^lLxN^ioL)~^'^NxL<^ subject to a^ct < P, (57) 

a 

to yield the error exponent in the AWGN case with correlated sensing noise: 

1 9"^ 

£awG-n{N) < — -^(Q:AWGN)^lixArR(aAWGN)~^lAfxLQ:AWGN) (58) 

where ckawgn provides to solution to (l57l) and the inequality in (l58l) is due to ([TT]) and the modified 
power constraint in (l56l) . To fully compute an upper bound on the right hand side of (|58l) . first, R(q;) is 
inverted and simplified. For the case of correlated noise, R(a) is given by 

R(a) = l^xLD(Q:)R^D(a)^lix7V + a^I^ 

iD(a)^D(Q:)Uxiv + cT^Iiv, (59) 

where A ^ B indicates that the matrix (A — B) is positive semi-definite. Using the Sherman-Morrison- 
Woodbury formula for matrix inversion. 



11 . . . . . -1 



R(o:) ^ -^In ^ItvxL 

0",, cr,. 



diag 



^LxNJ-NxL 



Amin I '3^1 1 / 



2 



Ilxtv^. (60) 



^2 



Invoking the Sherman-Morrison-Woodbury formula for matrix inversion once again. 



R(a)"' ^ ^Ijv - ^Mljvxjv, (61) 



where 

^ -\ An 
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due to the fact that ol^ol < P from (|56l ). 

By substituting (|6T| ) in (|57] ). the solution to (1571 ) is upper-bounded by the solution to 

maximize ol^Ilxlol subject to ol^ cx < P. (63) 

OL 

The value of ct that maximizes (|63l ) is the eigenvector corresponding to the maximum eigenvalue of I^xl. 
scaled to satisfy the constraint with equality. Substituting this in (l58l) . the bound in (|3TI) obtained, with 
the substitution, 7s = where 7s = 9'^/\ain and P-p < P/{pi6'^ + Amm)- 

(iii) Upper-bound (no sensing noise): With no sensing noise, = Olxl- The optimization problem to 
obtain the best error exponent is the same as in (|32l) . to yield (|36l ). 

Combining the modified AWGN upper-bound and the no sensing noise upper-bound in (|36l ). a joint 
upper-bound is obtained, which is identical to (|37l) . except for the substitution = Amin and % = 6'^/Amin- 
It follows that (|39| ) holds with z = 7c/ (pi7s + 1), to provide the proof. 
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